Testing the interaction model with cosmological data and gamma-ray bursts 
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\ We use the new gamma-ray bursts (GRBs) data, combined with the baryon acoustic oscilla- 

tion(BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release, the 
1 newly obtained A parameter at z — 0.6 from the WiggleZ Dark Energy Survey, the cosmic microwave 

^ 1 background (CMB) observations from the 7- Year Wilkinson Microwave Anisotropy Probe (WMAP7) 

o" 



O 

U 



results, and the type la supernovae (SNela) from Union2 set, to constrain a phenomenological 

model describing possible interactions between dark energy and dark matter, which was proposed 

to alleviate the coincidence problem of the standard ACDM model. By using the Markov Chain 

Monte Carlo (MCMC) method, we obtain the marginalized la constraints f2 m = 0.2886 ± 0.0135, 

r m = —0.0047 ± 0.0046, and wx = —1.0658 ± 0.0564. We also consider other combinations of these 

data for comparison. These results show that: (1) the energy of dark matter is slightly transfer- 

O f ring to that of dark energy; (2) even though the GRBs+BAO+CMB data present less stringent 

Q ■ constraints than SNe+BAO+CMB data do, the GRBs can help eliminate the degeneracies among 

S— I ' 

■ parameters. 
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INTRODUCTION 

A lot of astrophysical and cosmological observations have indicated that the universe is undergoing an accelerating 
expansion and great efforts have been made to understand the driving force behind the cosmic acceleration [l - 15 1 . This 
gives birth to the construction of a strange dark energy with negative pressure, which may contribute to interpret the 
present accelerated expansion. The most simple candidate for these uniformly distributed dark energy is considered 

\ 

to be in the form of vacuum energy density or cosmological constant (A). However, despite its simplicity, the simple 
cosmological constant is always entangled with the coincidence problem: Why the present matter density p m , which 
decreases with the expansion of our universe with a~ 3 , is comparable with the dark energy density pa, which does not 
change with the cosmic expansion of our universe? In order to relieve the coincidence problem, other dynamic dark 

3, phantom QQ, k-essence 20, 21 1, 



energy models were proposed in the past decades, including quintessence 
as well as quintom model 22 
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24]. However, the nature of dark energy is still unknown. The other presumption is 



naturally considered that energy is exchanged between dark energy and dark matter through interaction. 
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We assume that dark energy and dark matter exchange energy through interaction term Q, namely 



px + 3Hp x (l + w x ) = -Q, 
Pm + 3Hp m = Q, 



(1) 
(2) 



and the total energy conservation equation expresses as ptot + 3i? (ptot + Ptot) — 0, where p to t — Px + Pm- Because 
the format of interaction term still can not be deter mined f rom fundamental physics, many literature have extensively 
considered various forms of the interaction term Q 
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In this work, we consider the simplest case for convenience, i.e. 

Q = 3r m Hp 

III •■ 



38] 



(3) 



where r m is a dimensionless constant: r m = indicates that there is no interaction between dark energy and dark 
matter; the energy is transferred from dark matter to dark energy when r m < 0, and from dark energy to dark matter 
when r m > 0. We assume that the equation of state (EoS) of dark energy, wx = Px/ Px, is a constant in a spatially 
flat FRW metric. In this case, Eq. (JIJ reads 



px + SHp x {l + w x ) = -3r m Hp m . 
Combining Eq. ((2J and Eq. (jH), we can get 

Pm = Pm0(l + Zf^-^ - Prnoa' 3 ^-^ 



By using Eq. ([5]), we obtain 



px = Aa- 3(1+W ^ 



r mpmQ _3(i_ rTri ) 



(r m + w x ) 

where A is an integral constant. Inserting it into the Friedmann equation 



H 2 



K 2 (p m + Px)/3, 



(4) 



(5) 



(0) 



(7) 



where n 2 — 8ttG with G the gravitational constant, and this integral constant can determined with H(z = 0) = Hq 
The corresponding Hubble parameter is 



E\z) = (H/H Y 



(r m + w x ) 



(1 + *) 



+ ( 



(r m + wx) 



){l + z) 3 ^ 1+wx l 



(8) 



Recently the simple phenomenological interacting scenario has been constrained from several cosmological observa- 



tions 
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Also, the Gamma-ray bursts(GRBs) have been propos ed as distance indicators and regarded as a 
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55| . In contrast to supernovae, the high energy 



complementary cosmological probe of the universe at high redshifts 
photons in the gamma-ray band are nearly unaffected by dust extinction. Therefore, those observed high-redshift 
Gamma- Ray Bursts (GRBs) at 0.1 < z < 8.1 may constitutes are a complementary probe to fill the "desert" between 
the redshifts of SNIa and CMB. Nevertheless, due to the lack of a low-redshift GRBs at z < 0.1, there is a circularity 
problem in the direct use of GRBs j^]. Some statistical methods have been proposed to alleviate this problem, such 
as the scatter method 43 1, the luminosity distance method |43| and the Bayesian method {44]. Other methods trying 



to avoid the circularity problem have been proposed in Ref. [5g, |57| . Using the cosmology- independent calibration 



method proposed in Ref. 5 



59| have obtained 59 calibrated high-redshift GRBs called "Hymnium" GRBs sample 



out of 109 long GRBs with the well-known Amati relation. Therefore, it may be rewarding to test the coupling 
between dark sectors with this newly obtained 59 GRBs deprived of the circularity problem. 

With this aim, in this paper, we adopt the Markov Chain Monte Carlo(MCMC) technique to constrain one inter- 
action model from the latest observational data. To reduce the uncertainty and put tighter constraint on the value of 
the coupling, we combine the GRBs data with the joint observations such as the 557 Union2 SNela dataset Q|, the 



^results, an d the model independent 
601, the two BAO distance ratios at 



CMB observation from the Wilkinson Microwave Anisotropy Probe (WMAP7) 
new A parameter from the baryon acoustic oscillation (BAO) measurements 
z = 0.2 and z = 0.35 from the spectroscopic Sloan Digital Sky Survey (SDSS) data release7 (DR7) galaxy sample 
jsj . This paper is organized as follows: we introduce the observational data in Section [2] The numerical analysis 
results are discussed in Section [3J and the main conclusions are summarized in Section |U 



2. OBSERVATIONAL DATA 

The 59 calibrated Hymnium GRBs and the 557 Union2 SNela data sets are given in term of the distance modulus 
Theoretically, the distance modulus can be calculated as 

<" = 5lo STT~ + 25 = 51og 10J ff dL - Mo, (9) 
Mpc 



where /io = 5 \og 1Q [Ho/ (100km/ s/ Mpc)] +42 • 38, and the luminosity distance <1l can be calculated using cLl 

dz' 

W) ■ 



^TT^ So EU') ■ The X 2 value of the observed distance moduli can be calculated by 



where /i bs(zi) is the observed distance modulus for the SNe la or GRBs at redshift Zi with error a^; /J-(zi) is the 
theoretical value of the distance modulus calculated from Eq. ([9]). The nuisance parameter h is marginalized with a 



flat prior, after which we get 62 1 



^ a? (EtW<-lnlO/5) 2 



X 



(EtW<-lnlO/5) 2 nl /In 10 2V 



where = fi b s (zi) — 25 — 5 log 10 HqcL^. In this part the radiation component of the total density is neglected, because 
its contribution in low redshifts is negligible. 
_Fot the BAO observation, we use the new BAO A parameter at z = 0.6, with the measured value A = 0.452 ±0.018 
from the WiggleZ Dark Energy Survey. It can be calculated as the following: 



, _ /r -H D v (z = 0.6) _ yTi', 

J\ \j \L'XYi 



z = 0.6 0.6 
Now we can add our \ 2 obtained before with 



°' 6 'M^CmY" 



#(0.6) \n, 



(ii) 



A-0.452 N 2 



X 2 BAOa(P) = y 0~018 ) ■ (12) 

Notice that the BAO A parameter does not depend on the baryon density flbh 2 or the Hubble constant h. The 
radiation density which does depend on h is negligible to the Hubble parameter E(z). 
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In addition to the above A parameter, we also consider the BAO distance ration (d z ) at z = 0.2 and z = 0.35 from 
SDSS data release 7 (DR7) galaxy sample [fil|. The BAO distance ration can be expressed as 

j r a (z d ) 



where the effective distance is given by 11 1 

Dv(z) 

and the drag redshift Zd is fitted as [f^ 



D v {z)' 



d\{z) z i 1/3 



{l + zYH{z) 



(13) 



(14) 



1291 (n m /i 2 ) 0251 , n 9sh , 

6 a = 0.313(fi m /i 2 )- 0419 [l + 0.607(O ro /i 2 ) 0674 ], 

b 2 = 0.238(fl m /i 2 ) ' 223 , (16) 

The comoving sound horizon is 

[°° c s (z)dz 

rs{z) = J z -eW' (17) 



where the sound speed c s (z) = 1/ x/3[i + i? b /(l + i)], and i? b = 3ft b /i 2 /(4 x 2 - 4 69 x 10" 5 ). The x 2 value of BAO 
observation can be expressed as [6l| 

Xbao = APb AO Cbao _1 AP B ao, (18) 

where APbao = Pth — Pobs, Pobs is the observed distance ration. Cbao 1 is the corresponding inverse covariance 
matrix. 

For the CMB observation, we use the derived dataset from the WMAP7 measurement, including the acoustic scale 
(Z a ) ,the shift parameter R, and the redshift of recombination z* ij]. The acoustic scale can be expressed as 

«a=f 7 — v , (iyj 

where r s (z*) = Hq^ 1 c s (z)/E(z)dz is the comoving sound horizon at photon-decoupling epoch. The shift param- 
eters can be expressed as 

R(z.) = ^SLinn ( ^\Sh\ t ■ (20) 



The decoupling redshift is fitted by 



G4 



VW\ \ 'Jo E(z) 



z« = 1048[1 + 0.00124(fi b /i 2 )-°- 738 ][l + gi(n m h 2 ) 32 ], (21) 

_ 0.0783(r! fc fe 2 )- - 238 _ 0.560 
91 ~ l + 39.5(fW) 0763 ' 92 ~ 1 + 21.1 (n b h 2 ) 181 ' ( ' 

The x 2 for CMB can be expressed then as 

Xcmb = APJ mb Ccmb _1 AP C mb, (23) 

where APqmb = Pth — Pobs, Pobs is the observed dataset, and the Ccmb 1 is corresponding inverse covariance 
matrix. 
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Data 


wx 




O m 


H / 100 


X 2 /dof 


all 


-1.0658 ± 0.0564(l<r) 


-0.0047 ± 0.0046(1<t) 


0.2886 ± 0.0135(1ct) 


0.7369 ± 0.0391(1ct) 


0.9140 


all-GRBs 


-1.0657 ± 0.0562(ler) 


-0.0047 ± 0.0046(lcr) 


0.2886 ± 0.0135(1ct) 


0.7367 ± 0.0389(ler) 


0.9686 


all-SNe 


-1.2757 ± 0.2358(1<t) 


-0.0062 ± 0.0045(1ct) 


0.2649 ± 0.0288(lcr) 


0.7984 ± 0.0803(1<t) 


0.3580 



TABLE I: The marginalized la errors of the parameters wx, r m , f2 m , and Ho for the phenomenological scenario, as well as 
X 2 /dof, obtained from the combinations of the data sets GRBs+SNe+BAO+CMB (all), SNe+BAO+CMB (all-GRBs), and 
GRBs+BAO+CMB (all-SNe). 

3. CONSTRAINT ON THE PHENOMENOLOGICAL INTERACTING SCENARIO 



The model parameters are determined by applying the minimum likelihood method of \ 2 fit- Basically, the model 
parameters are determined by minimizing 

X 2 = XGRBs + XSNe + XcMB + XbAO + XBAOa- ( 24 ) 



We apply the Monte Carlo Markov Chain (MCMC) method 65[ with 8 chains and obtain the marginalized la 
constriants: Q, m = 0.2886 ± 0.0135, r m = -0.0047 ± 0.0046, and w x = -1.0658 ± 0.0564. We show the marginalized 
probability distribution of each parameter and the marginalized 2D confidence contours of parameters in Figure[T] For 
comparison, fitting results from the combinations of SNe+BAO+CMB and GRBs+BAO+CMB are shown in Figure 
[2 Because CMB and BAO data tightly constrain the matter density f2 m , and therefore help improve the constraints 
on dark energy property [15I ] , so these data are taken as priors in our treatment and are combined with other data to 
test the constraining power of GRBs and SNe data. It is shown that the plots in Figure Q] and the left panel of Figure 
[5] are almost the same, suggesting that the current GRBs data are consistent with other observations, although they 
contribute little to the existing constraints. Comparing the constraints from SNe data with those from GRBs data 
in Figure [2l we see that the degeneracies between Q m and wx , Q m and r m are different although the current GRBs 
data give larger errors on wx- Therefore, GRBs data have the potential to help constrain the model parameters wx 
and r m by combining with SNe data. The best-fit values of the parameters along with their la uncertainties from all 
three different combinations mentioned above, are explicitly presented in Table [I] In general, for the three different 
joint data sets, at the la confidence region, the energy is seen transferred from dark matter to dark energy, and the 
concordance ACDM {wx = — 1, r m = 0) model can not be excluded. 

From Figure [TJ we see that the coupling parameter r m is correlated with all other parameters f2 m , wx and 
Hq. However, r m = is within the la confidence region. To see the point clearly, we fix r m =0 and obtain 
the marginalized la uncertainties of the other model parameters from the combination of all observational data 
(SNe+CMB+BAO+GRBs): Q m = 0.2863 + 0.0130, w x = -1.0436 + 0.0509, and # /100 = 0.6994+ 0.0118. The 
constraints on the model without interaction (r m = 0) are consistent with those with interaction at la level, and Wx 
shifts toward —1 when r m = so that it is consistent with ACDM model at la level. 



4. CONCLUSIONS 



In this paper we constrain an interacting dark energy model using the new GRBs, the Union2 SNe, CMB, and 
BAO data sets. By adopting the MCMC approach we obtained the marginalized la errors of each parameter: 
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FIG. 1: The marginalized probability distributions, and the marginalized la and 2a confidence contours of the parame- 
ters wx, r m , flm, and Ho in the phenomenological interacting scenario, from the combinations of all observational data 
SNe+BAO+CMB+GRBs. 




FIG. 2: The same as in Figure [T] but the constraints in the left panel are from the data combinations SNe+BAO+CMB 
(all-GRBs) and the constraints in the right panel are from the data combinations GRBs+BAO+CMB (all-SNe). 
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Q m = 0.2886 ±0.0135, r m = -0.0047 ± 0.0046, and w x = -1.0658 ± 0.0564. 

In order to test the constraining power of the observational GRBs data, we compared the results from the combina- 
tions of GRBs+BAO+CMB with those from SNe+BAO+CMB. We find that the constraints from the current GRBs 
data are less stringent than those from SNe data. The error bars of wx, and Hq from GRBs+BAO+CMB are 
roughly twice of those from SNe+BAO+CMB, although the error bar of r m is roughly the same as shown in Table U 
and Figure O However, the directions of degeneracies between f2 m and W X > Q m and r m are different, so GRBs data 
have the potential to help tighten the constraint on the parameters O m , wx and r m if the measurement precision 
of the data is improved in the future. By fitting the model to all the observational data combined, we find that the 
energy slightly transfers from dark matter to dark energy at la region. We also note that the coupling parameter r m 
is correlated with all other model parameters fi m , wx and Hq, and r m — is within the la confidence region. The 
constraints on the model without interaction (r m = 0) are consistent with those with interaction at la level, and the 
value of wx shifts up a little when r m =0. In conclusion, the concordance ACDM model still remains a good fit to 
the observational data. 
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